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Abstract 

We obtain analytical expressions for the generating function £d(A) 
of the sum of ci-squares arithmetical function r<j(n) where A is a free 
parameter. The original d-dimensional infinite sum is reduced to a 
formula containing a single finite sum over a convergent series. We 
compare the formulas to numerical computations and show that the 
percentage difference is negligible at small A for various values of d. 
£d(A) divides naturally into two terms and we show that one term has 
a direct physical application to the <i-dimensional Casimir energy of 
massless scalar fields in cubic cavities. 
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I. Introduction 



The representation of integers as a sum of d squares is one of the classical 
problems in Number Theory pQ . The number of representations of a positive 
integer n by d squares, allowing zeros and distinguishing signs and order, is 
denoted by the arithmetical function r^(n) where it is common to write r(n) 
for r2(n). Great ingenuity and effort has gone into obtaining explicit formulas 
for these functions. In particular, a systematic treatment for even values of 
d < 8 existed by 1940. However, a systematic treatment for odd values 
with d < 7 was only obtained in 2002 by Shimura [2]. Explicit formulas for 
any value of d do not presently exist. Fortunately, explicit knowledge of the 
functions is not required for the calculation of generating functions. In this 
paper we are interested in a particular generating function of the arithmetical 
function rd(n) because of its connection with a vacuum boundary effect in 
quantum field theory called the Casimir energy. The most-well known cases 
of infinite series of r^(n) are the Dirichlet and power series involving r(n). 
The Dirichlet series is given by PQ, 

OO , s 

Y = 4 COO L(s) for s > 1 where L(s) = V - 3~ s + 5~ s - . . . (1) 
^— ' n s 

n=l 

while the power series is given by pQ, 

oo 3 5 

Y r (n) x n = 4 ( — X — + -..). (2) 

n=l 

An important result is that the average order of r(n) is simply tt: 

N 



l^r(n) = 7r. (3) 



hm 

N^oo N 

n=0 

There is a clear geometrical interpretation of the above result. The sum up 
to N is equivalent to counting the lattice points inside and on a circle given 
by x 2 + y 2 = r 2 where the radius r = y/~N. If each lattice point is assigned 
an area of one unit, the total area approaches the area nr 2 and hence tt N 
as r approaches infinity. Good estimates of the correction to it N for finite 
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N have been obtained as well as the distribution of the lattice points on the 
circle 

In this paper, we are interested in results where the dimension d is not 
fixed a-priori. Specifically, we obtain accurate formulas for the infinite series 
of the arithmetical function r^(n) multiplied by an exponential damping term 
raised to the square root of n. The generating function is: 

oo 

aA)E^r4n)e-^ (4) 

n=0 

where the free parameter A is a positive real number and yields a convergent 
series. £d(A) is a function of both the parameter A and the dimension d. Note 
that it deviates from usual power series like (J2J) because the exponential is 
raised to the square root of n instead of n. The \Jn is chosen for physical 
reasons. The frequency of normal modes of a massless scalar field in a d- 
dimensional cubic cavity is proportional to \/n (where n = n\ + n% + ■ + n 2 d 
and rii is the mode in the ith direction). These can then be summed in a 
particular fashion to obtain the <i-dimensional Casimir energy. Sums over 
the \fn can be obtained from (J3J) by taking the derivative with respect to 
—A. More generally 

00 ( r\\ v 

J2r d (n)n^e- x ^=(-iy - £ d (A) (5) 

n=0 ^ ' 

where p is a non- negative integer. Hence, with p — 1, we have a sum over 
the \fn and a possible connection with the Casimir energy in cubic cavities. 

The formula we obtain for £d(A) divides naturally into two terms: an inte- 
gral term id (A) and a second term which we label Cd(X). id(A) is equivalent 
to replacing the multiple sums in £d(A) by multiple integrals. For small A, 
irf(A) can be viewed as a first approximation to £d(A) with C^(A) a small 
correction. However, at large A, Cd(A) actually dominates over id(A). We 
compare our formulas to direct numerical computations for values of d rang- 
ing from 2 to 5 and for A ranging from 0.1 to 10. The difference (error) is due 
to a remainder term we neglect in the derivation of our formulas. For small 
A, the agreement is excellent. For A < 1 we obtain less than 0.02% error and 
for A < 5 we obtain less than 6.5% error. The percentage difference increases 
as A and d increases. At A = 10 the error is 10% at d = 2 and 32% at d — 5. 

The remarkable thing is that the formula for £d(A) is already in a form 
that contains the necessary pieces for a clear formulation of the Casimir 
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energy. In the end, all the information we need is contained in Cd(A). We 
show how one can obtain the d- dimensional Casimir energy of massless scalar 
fields in cubic cavities by taking the appropriate derivatives and limits of the 
parameter A in Cd(A). 



II. Formulas Obtained 

In this section we are interested in obtaining analytical expressions for the 
generating function 

oo 

a^E r ^ e " iVH ' ( 6 ) 

n=l 

The above definition is equivalent to 

oo oo 

&(A) = -1+ ■•• E e- . (7) 

nj=— oo ni=— oo 

Define the following short-hand form for a j-dimensional multiple sum: 

oo oo 

(EV = E ■•• E e - A v / "i+«i+-+« J 2 . ( 8 ) 

nj=— oo m=— oo 

A sum from minus to plus infinity that excludes zero is labeled . It is 
convenient to construct the following (j + l)-dimensional sum: 

00 oo 00 

= E' E ••• E e-v» a ^^-^. ( 9) 

?i=— 00 rij=— 00 m=— 00 

Note that the last sum in Aj(A), the sum over the variable labeled n, does 
not contain a zero. The d- dimensional sum in (JJJ) can now be expanded in 



4 



the following fashion: 



^+£+£'£+£(£) 2 +- + £(£ 
= i+££'(£) J 



(10) 



i=o 

d-1 

= l + J>i(A) 

i=o 

Substituting (JTUD into (J7|) yields 



d-1 



6^) = ^ A, (A). (11) 

j=0 

The goal now is to find analytical expressions for Ay (A). To this end we 
apply the Euler-Maclaurin formula that converts sums to integrals. The 
Euler-Maclaurin formula is [3]: 

00 poo 1 1 1 

£/(")=/ /Wdx--/(0) + ^ — B 2p /M(0)+fi 5 (12) 

n=l ^° p=l ^ V'' 

where /( 2p_1 )(0) are odd derivatives evaluated at zero and q is a positive 
integer. R q is the remainder term given by [3] 



1 pi 00 

^ = -Tom / E ^ ( x + ") rfa; 

^ ! Jo u=0 



(13) 



where B2 q (x) are Bernoulli functions. We will neglect the remainder term 
R q in deriving formulas for £<f(A). The error introduced by omitting R q 
is calculated in the section where we compare the formulas to numerical 
computations. 

We want to make use of the Euler-Maclaurin formula to determine Ay (A). 
The function / being summed is the exponential function in ©. For this 
function, / 2p_1 (0) is zero for the first j sums in (jOJ) (all sums except the last 
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one). A proof of this is given in [7j. When / 2p_1 (0) = 0, the sum from p = 1 
to q in(|12j) is zero independent of q. This implies that R q given by (|13J1 must 
be independent of q for the case of our exponential function. We prove this 
explicitly in Appendix A. Note that if R q is nonzero, it has the same nonzero 
value for all q even though l/(2g)! approaches zero as q — > oo. 

With / 2p-1 (0) = for the first j sums we neglect R q and the Euler- 
Maclaurin formula (J 12)) for the first j sums reduces to 



°° PC 

£/(») = / 

n=l J 



/(x)dx-i/(0). (14) 



The function / in @ has the property fijii) = /(— n»). The sum over a 
given nj can therefore be written as 

oo oo 

f(jn) = 2 /("*) + /(0) 

n;=— oo rii=l 

= 2N°°/(x)^-i/(0)j +/(0) (15) 



/(x) 

where we used (JHJ). From (|T3j) we see that each sum in @, except the last 
one, can be replaced by an integral. This yields 

•V A: E(E) ' Y. f 

00 . poo , 

= J2 e- x V n2 +^+---^dx 1 dx 2 ...dx j 

n=-oo 

°° ^oo 

= 2- 7+1 V / e- A V n2 +*l+*i+~+*Ux 1 dx 2 ... d Xj 

n=l ^ 



(16) 



The integral can be expressed in terms of a modified Bessel function 

V"--' ■■%!.,-,... d,- 



C~ 



(17) 

- 2 v^ (%i(An) (^ 
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The modified Bessel function Kj-i (A n) has the following identity j3j: 

K^Xn)(lf^(-l)^(jLfK 0{ X n) . (18) 
Substituting flUD and JT7J) into (HHJ) yields 



A,(A)=2^^(-1)^AM )^ J> (An). (19) 



1 oo 



d\\\d\ 



n=l 



The infinite sum over the modified Bessel function K (Xn) has the following 
series expansion [3]: 



oo oo / -, i \ 

VK (An) = -{C + ln(A/47r)} + ^- + 7ry ( - J 

^ 1 ; 2 1 1 7 ;/ 2A ^ Vv / A 2 + 4m 2 7r 2 2 run J 

n — 1 m— 1 \ v ■ y 



1 VVA 2 + 4m 2 7r 2 2m7r 

n=l m=l x v 

Substituting flUD into (HH yields 



(20) 



A,-(A) = ~V vr 2 ^ r(i±i) + -4t 2 j+1 tt* r(^) 
n A^ v 2 A J+ v 2 (21) 

+ A2 i +2 r(if) 7T* Xi(A) 

where 

oo 1 

^) s £ , u , a a ^ • ( 22 ) 
m=1 (A" 1 + Am z it 2 ) 2 

Substituting (0} into (HH) yields: 
£d(A) = E -| ^ r(i±i) + |^ vri T(2± 2 ) + A 2^ 2 r(^± 2 ) vri Xj (X) . 

3=0 

(23) 

Three terms are being summed. Note that if the first term is written as — o 3 -, 
then the second term is Oj+i. The sum over j of the first two terms is then a 
telescopic sum equal to — a + a d = — 1 + |^ 7r~ T(^-). Substituting this 
result into (J23|l yields 

&(A) = S ^ r (¥) - 1 + A E 2 J+2 T(^) tt* Xj(A) . (24) 

3=0 



Equation (|24j) is our final formula for £d(A). The original <i-dimensional 
infinite sum ((7j) has been replaced by a finite sum over the convergent series 
Xj(A) given by (|22j). It is convenient to express (}2"4")) in the form 

6i(A) =i d (A)-l+C d (A) (25) 

where 

W) = ^ vr^ r(^) and C d (A) = A ]T 2^' +2 T(^) tt* Xi (A) . (26) 

5=0 

There is a clear interpretation to (J23j). id(A) is identified as the 'integral con- 
tribution' to £<j(A). It is the multiple integral of e~ Ar over the <i-dimensional 
volume where r = (x\ + x\ + • • • + x 2 d ) l l 2 . This is shown in Appendix B. 
The expression id(A) — 1 appearing in (f2"5|) is equivalent to (jZJ) with the <i 
sums replaced by d integrals. This can be viewed as a first approximation to 
£d(A). Cd(A) is the interesting part for mathematical and physical reasons. 
Mathematically, C^(A) is the non-trivial 'correction' to the approximation of 
replacing sums by integrals in (j7J). Physically, C^{X) is the relevant part for 
the calculation of the Casimir energy as we will demonstrate later. Cd(X) is 
not necessarily smaller than /^(A). Their limits as A —>■ and A — > oo are: 

C<i(A) — > and id(A) — > oo as A — > 

and (27) 
Cd(A) — * 1 and /^(A) — > as A — >• oo . 

It is clear that id (A) dominates at sufficiently small A and Cd(A) dominates 
at sufficiently large A. How small is sufficiently small? For a given dimension, 
this depends on A. Values are quoted in Table 1 for A ranging from 0.1 to 
10 and d ranging from 2 to 5. In this sample, i^(A dominates for A = 1 and 
below and Cd(X) dominates for A = 5 and greater. It is therefore only valid 
to view Cd{\) as a correction term for A below a certain value. 

Formula (|24j) together with ()22|) for Xj(A) are our final formulas needed 
to determine £d(A). We can go one step further if 0<A<27r. Then XjW can 
be expanded as a power series in A involving Riemann zeta functions. The 
result is a power series for £d(A) which converges rapidly when A << 2n. 

Power Series Expansion for Xj(A) 

For the case < A < 2n, Xj(A) can be expressed as a power series in A with 
Riemann zeta functions as coefficients. We begin by rewriting (122)) in a form 
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Table 1: Values of Crf(A), /d(A) and their Ratio 





Lambda 


'Correction Term': C d (X) 


Integral Term: I d (X) 


Ratio: C d (X) / l d (X) 


D=2 


0.1 


0.022751 


628.318531 


3.62 E-5 


1 


0.222914 


6.283185 


3.55 E-2 


5 


0.778569 


0.251327 


3.10 


10 


0.937333 


0.0628324 


14.9 




D=3 


0.1 


0.026241 


25132.741228 


1.04 E-6 


1 


0.256219 


25.132741 


1.02 E-2 


5 


0.84941 


0.201061 


4.22 


10 


0.9751 


0.025128 


38.8 




D=4 


0.1 


0.0287472 


1184352.5281307 


2.43 E-8 


1 


0.2798216 


118.4352528 


2.36 E-3 


5 


0.88810 


0.189494 


4.69 


10 


0.9884 


0.011882 


83.2 




D=5 


0.1 


0.030834 


63165468.166972 


4.88 E-9 


1 


0.299219 


631.654682 


4.74 E-4 


5 


0.912840 


0.202130 


4.52 


10 


0.9940 


6.3299 E-3 


157 



suitable for a binomial expansion: 



*(A) = E 



m- 



\ (A 2 + 4 m 2 7r" i ■-• 



1 y Wi A 2 

(2ttV+ 2 ^ mi+A + 4m 2 7T 2 

v ' m=l 



(28) 



We can obtain a binomial expansion for ^1 + 4TT ^2 7r 2 ^ 2 if A 2 /4m 2 7r 2 < 1 

for all m. This is true if < A < 27r (the condition A > is there because 
A was assumed from the start to be a positive real number). The binomial 
expansion yields 



\2 . -U+2) °° , -(3+2) X 2 

4m / 7r" ! / -^^V n /\4m / 7r z 

n=0 



-l) n A 



2» 



oo 

E n ! 4^ y^n ^ 

n=0 i=l 



1 ™ 

srII(«' + J/2). 



(29) 
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(30) 



Substituting (JH into (J2ED yields 

I 00 (— l) n A 2n n 00 1 

= (2vrV+ 2 11^ + ^/ 2 ) E m i+2+2n 

^ ' n=0 i=l m=l 

v 1 n=0 ' i=l 

Expanding ()30|) term by term yields 

MX) = j^jr 2 {cU + 2) - ^ C(j + 4) (1 + j/2) 

+ 16^2! C(J + 6) (2 + j/2) (1 + j/2) + • • j • 

(31) 

Sustituting (HUD into (J2D yields 



2^ 

+ A E ^§{C0' + 2 ) - ^ C(7 + 4) (1 + J/2) (32) 

j=0 ^ 



+ CO' + 6) (2 + 3/2) (1 + j/2) + . . 

The above power series for £d(A) is valid when < A < 2n. One of its 
advantages is that it converges rapidly when A/27T << 1. 



II. Numerical Algorithm and Comparison to Formulas 

In the process of deriving a formula for £d(A) we omitted the remainder term 
R q . The formulas will therefore contain some errors. Recall that the exact 
expression for £d(A) is 

00 

aA) = $> d (n)e- A ^ 

00 00 (33) 

X ^ X ^ _ \ . /n2j_„2 j i_„2 



7id=— 00 ni 



^2 e ~ x V n i +n 2 + - +n d 
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We developed an algorithm that computes (|33|) in a time efficient way by 
making use of standard permutation formulas to avoid using negative inte- 
gers and all possible orderings. The algorithm is also constructed with the 
flexibility of allowing the dimension d to be a variable and user- input. In con- 
trast, if the numbers were generated using d- nested for- loops, then one would 
need to change by hand the number of for-loops for different dimensions. 

We then calculate £d(A) using the formulas we derived: (J23J) together with 
(|22|) (it is convenient to use the power series for cases when A is small i.e. 
A«27r). We then determine the absolute and percentage difference between 
the formula and algorithm. We performed calculations for the following cases: 
A = 0.1, 1, 5 and 10 and dimensions d = 2, 3, 4 and 5. The results are quoted 
in table 2. The numbers are quoted to the precision reached with the last 
digit rounded. 

For small A, the agreement is excellent. For A < 1 we obtain less than 
0.02% error and for A < 5 we obtain less than 6.5% error. It is important 
to note that the agreement at values like A = 5 is highly dependent on 
the contribution made by C^(A). As previously mentioned, Crf(A) makes 
a dominant contribution already at A = 5 (see table 1). The percentage 
difference increases as A and d increase. For the highest value of A in table 2, 
A = 10, the error is 10% at d = 2 and 32% at d = 5. Note that the absolute 
difference (the remainder term) changes only slightly with A in comparison 
to the percentage difference. 



III. Casimir Energy and the Generating Function of 

rd(n) 

Formula (j2*%Jl for the generating function £<z(A) of rd(n) has a direct rela- 
tionship to the Casimir energy of a massless scalar field confined to a d- 
dimensional box of side L. In fact, the natural division of £d(A) in two terms, 
Jrf(A) and Crf(A), is 'tailor made' for extracting the Casimir energy. Sim- 
ply put, the expressions for the regularized vacuum energy with and without 
boundaries are expressed in terms of £d(A) and /^(A) repectively. The Casimir 
energy is their difference and is written in terms of Crf(A). 

Consider a massless scalar field <p(x) confined to a box of side L with Neu- 
mann boundary conditions where the derivative of the scalar field vanishes 
at the boundaries i.e. di<ft(x) = when x,- t = or L. We now quickly go over 
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Table 2: Formula versus Numerical 





Lambda 


Formula 


Numerical 


Absolute 
Difference 
(Remainder Term) 


Percentage 
Difference 


D=2 


0.1 


627.341282 


627.341408 


1.26 E-4 


2.0 E-5 


1 


5.506099 


5.507242 


1.10 E-3 


2.0 E-4 


5 


2.9896 E-2 


3.0648 E-2 


7.52 E-4 


2.5 


10 


1.6539 E-4 


1.8449 E-4 


1.91 E-5 


10.4 




D=3 


0.1 


25131.767469 


25131.767881 


4.12 E-4 


1.6 E-6 


1 


24.3889599 


24.3926827 


3.72 E-3 


1.5 E-2 


5 


5.0471 E-2 


5.2752 E-2 


2.28 E-3 


4.3 


10 


2.282 E-4 


2.813 E-4 


5.31 E-5 


18.9 




D=4 


0.1 


1184351.5568779 


1184351.5566139 


2.64 E-4 


2.3 E-8 


1 


117.7150744 


117.7232799 


8.21 E-3 


7.0 E-3 


5 


7.7594 E-2 


8.2251 E-2 


4.66 E-3 


5.7 


10 


2.822 E-4 


3.8153 E-4 


9.93 E-5 


26.0 




D=5 


0.1 


63165467.197806 


63165465.850354 


1.35 


2.1 E-6 


1 


630.953901 


630.969211 


1.53 E-2 


2.4 E-3 


5 


1.14970 E-l 


1.22960 E-l 


7.99 E-3 


6.5 


10 


3.299 E-4 


4.855 E-4 


1.56 E-4 


32.1 



the standard steps leading to to the expression for the regularized vacuum 
energy. The vacuum energy is the sum over all modes of ~ u (we work in units 
where h = c = 1). For a massless scalar field confined to a ci-dimensional box 
of side L with Neumann boundary conditions, the frequency uj is given by 

u=l{n\ + nl + . + nW (34) 

where rii > 0. The vacuum energy is 

oo oo 

^ = ^E---5>? + ^ + - + ^) 1/2 (35) 

n^=0 ni=0 

which is formally infinite. We regularize it with an exponential term so that 
the regularized vacuum energy is 

oo oo 

H d (X) = ^E-E^ + n ' + - + ^) V2 e-^+^+-+^ 1/2 . (36) 

71^=0 ni=0 

We will later take the limit as A — > 0. For now, notice the resemblance 
between H d (X) and —d\£dW- Of course, the two are not equal since the 
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sum in £d(A) goes from — oo to oo. However, we can express the d sums from 
to oo found in Gd{X) in terms of sums from — oo to oo. The function we 
are summing in (J3T?)) has the property f(n) = /(— n) where n can be rii,n 2 , 
etc. We therefore have the relation /(n) = h J2°?oo f( n ) + I which 
can be expressed as an operator — > | (X^oo + -0 • Applying the operator 
d times yields 

oo, ^ oo,^d 7 oo 



P =i 



where the case p = yields zero i.e. /(ni = 0, n 2 = 0, = 0) = 0. We 
can therefore express (j3T?j) as 



Equation (j38j) shows already a connection between the generating function 
of rd(n) and the regularized vacuum energy. We now want to calculate the 
Casimir energy. The Casimir energy is the difference between the vacuum 
energy with boundaries and the vacuum energy without boundaries. Hd(\) 
represents the energy with boundaries. Without boundaries the fields are 
continuous: there are no discrete modes. This is equivalent to replacing 
sums by integrals. We have already seen that £ P (A) — > I P (X) — 1 when sums 
are replaced by integrals and hence, d\£, p (\) — > d\I p (\). So the regularized 
energy without boundaries, Gd(X), is 

G -W = -l2^E(p) 5 A/p(A). (39) 
p=i y 

To return to the original vacuum energy we need to take the limit as A — > 0. 
In this limit, both Gd(X) and Hd(X) are divergent. However, their difference 
in this limit is finite and equal to the Casimir energy Ed, 

7T 1 v^/(i N 



E d S \im(H d (\) - G d (\)) = -!—•£( ) Mm 9xfe(A) - I r (X» 

(40) 

p=l 



7T 1 \— K/d^ 

L2 d + l 
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where was used. As we mentioned previously, C P (A) is the only A — 
dependent term that finally enters into the calculation of the Casimir energy. 
From the definition ()26|) of Cd{\) we obtain 

p-i 

lim d x C p (\) = $>+ 2 nf)^ lim Y,(A) 

A — *U A — *U 

= £r(ffi)^ C (j + 2) 

3=0 

where we used the power series fJ3T]) for XjW to obtain 

Substituting (J4T]) into (j4Ti|) yields the Neumann Casimir energy for massless 
scalar fields in a box of d dimensions: 

^4^£E(>¥)^C0 + 2 ). (43) 

p=l j=0 F 

Note that equation ()43|) shows clearly that the Neumann Casimir energy 
is negative since the terms inside the sum are all positive. Values calculated 
using (J4~3j) are in agreement with previous values calculated for the Neu- 
mann energy using different techniques (HI El QH] ■ Similar expressions to (J43|) 
can be obtained for other boundary conditions such as Dirichlet boundary 
conditions. 
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A Appendix 

In this appendix we show that the remainder term R q is independent of the 
positive integers q for the exponential function 



where C and A are positive real numbers. We want to show that R q = R q +\ 
where 



The function f(x) has the following properties: 

• odd derivatives evaluated at zero are zero i.e. / 2p-1 (0) = 0. This is 
proven in [7j. 

• the function and its derivatives are zero asymptotically 
i.e. lim^oo f(x) = and lim^-a, f p (x) = 0. 

The Bernoulli functions have the following properties: 



f(x) = e 



XVx' 2 +C 



(44) 




(45) 



• B'Jx) = —B n (x) =n5„_i(i). 

ax 

• B n (l) = (-1)" B n (0). 

• B 2n+1 (0) = except for ^(0) = -1/2. 



We use the above properties to obtain the following equality: 



B 2q (x) {x + v) = -f^ (x + v) 




(46) 



The integral in ()45j) is then given by 




(47) 
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where we used the property B 2q+ i(l) = B 2g+ i(0) = 0. Substituting (|47|) into 
(gnj yields 



R n 



(2q 



-— / B 2q+1 (x) f q+1 (x + v) dx . 



(48) 



We now repeat the process one more time. We obtain the equality 
B 2q+1 (x) (* + „) = f 2 ^ (x + u) 

= ±rB 2s±2 {x)_ 

dx\ 2q + 2 J K 
_ B 2q+2 {x) 2q+2 

2q + 2 J K ' 

The integral in ()48|) is then given by 

f B 2q+1 (x)f 2 ^ (x + u)dx = f 2 ^ (x + v) 

Jo 2q + 2 



(49) 



(50) 



2g + 2 J 



B w (x)f« +2 (x + u)dx 



The first term yields, 



2g + 2 



B^aCl^Cl + i/J-S^aCO)/ 2 ^ 1 



(51) 



= B 2q+2 (0) (f« +1 (1 + ,/) - (v)) 
where we used the property S 2g+2 (1) = -829+2(0). Summing over z/ yields, 



£>w°) (/ 29+1 a + - / 29+1 w) 



y=0 



£ 2g+2 (0) ( Km / 2 * +1 (1 + V ) - f« +1 (0)) = . 



(52) 



where we used the properties of f(x) previously listed. Therefore, summing 
the first term over v in (joTIj) yields zero. Substituting (J5Uj) into (J4*%|) yields 



(2g + 



-1^ / i?2 g+ 2(x)^^ +2 (x + Z,)rfx 



(53) 
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which is the desired result. Since R q = R q+ i, the remainder term is indepen- 
dent of q. 



B Appendix 

In this appendix we show that the term, 2 d r(^±i)/A d , which appears 
in the formula (J24j) . is a d- dimensional volume integral obtained by replacing 
the sums by integrals in (JJJ). Our goal is to show that 

/CO cyd 

We switch from cartesian to spherical coordinates (r, 0, 9i, 9 2 , . . . , ^-2) where 
r = (xf + ^2 + • • ■ + 3 ; rf) 1 ^ 2 - The d- dimensional volume element dV in spherical 
coordinates is 

dV = r d ~ l dr dcf) sin 9 1 d9 l sin 2 9 2 d9 2 . . . sm d ~ 2 d _ 2 d9 d _ 2 

d ~ 2 (55) 
= r d - 1 drd(j)Y[sm k 9 k d9 k 

k=l 

where 2 7r > > 0, 7r > 0j > and r runs from to 00. Substituting (jSTjj) into 
(J53) yields 

I d (X)=2TT e _Ar r*- 1 dr / [}sin fe fc d0 fc (56) 

where the 2 7r from the integral over was factored out. The integral over 
any 9 yields 

/ sin k 9 h d9 h = ^-^ (57) 
Jo 1 (—J 

and therefore 

d-2 d-2 



The integral over r yields 



00 



e -Ar r <i-l dr= EM. (59) 
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Substituting (J3HD and ® into (JUJ) yields 



I d (X) = 2 



?ri r(d) 
r(d/2) A d 



(60) 



The above expression can be simplified if we use the identity 



r(d) 1 2 d r(^±i) 
r(d/2) ~ 2 7^ ' 

Substituting (JBTj) into (joTIj) yields our desired result 



(61) 



2 d 

>7 



vr^ r(^±i). 



(62) 
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